Breakdown of perturbative weak coupling approaches for the biomolecular energy 

transfer 
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We show that the biomolecular exciton dynamics under the influence of slow polarization fluctu- 
ations in the solvent cannot be described by approaches which are perturbative in the system-bath 
coupling. For this, we compare results for the decoherence rate of the exciton dynamics of a resumed 
perturbation theory with numerically exact real-time path-integral results. We find up to one order 
in magnitude difference in the decoherence rate for realistically slow solvent environments even in 
the weak coupling regime, while both results coincide for fast environmental noise. This shows 
cxplicitely the nonperturbative influence of the bioenvironmental fluctuations and might render 
current perturbative approaches to biomolecular exciton transport questionable. 



I. INTRODUCTION 

Photosynthesis in biomolecular complexes is a funda- 
mental but at the same time rather complex series of 
processes which converts photonic energy of sunlight into 
chemical energy which is stored by an organism for later 
use to drive cellular processes 0,0. At the be gmmng 
of this chain, a photon is absorbed by one of the light- 
harvesting antenna structures (pigments) and an exciton 
is created. Its excitation energy is subsequently trans- 
ferred via several intermediate pigments within the an- 
tenna system to the reaction center where it is converted 
into chemical energy. The transfer of this excitation en- 
ergy within the antenna structure occurs by migration of 
electronic excited states from one molecular part to an- 
other and is a purely physical process which astonishingly 
occurs with almost unity quantum yield. 

The theoretical formulation for energy transfer which is 
nowadays widely applied to photosynthetic systems was 
provided by Theodor Forster in the 1940s 0]. It holds for 
the case when the pigments arc spatially rather distant 
and thus weakly coupled by incoherent (diffusive) exciton 
flow and also when they are very close and thus strongly 
coupled by coherent exciton motion The latter is 
often referred to as exciton coupling and the resulting 
exciton-coupled dimer can effectively be viewed as a su- 
permolecule with dclocalized electronic transitions where 
the molecular wave functions are given as superposition 
of the wave functions of the two monomeric pigments. 
Despite the somewhat different situations, it has been 
shown that the two cases are just two sides of the same 
coin 0, H, H| . The Forster energy transfer mechanism is 
a nonradiative resonance transfer process based on the 
Coulomb dipole-dipole interaction of the excitons in the 
pigments. 

The photosynthetic excitons are by no means isolated 
quantum objects, but are embedded in host protein struc- 
ture and in the surrounding polar solvent. It has been 
shown Q that decoherence due to the coupling of the 
electric dipole moment of the exciton to the fluctuating 
electric dipole moments of the individual molecules of 



the solvent is not negligible. Decoherence and energy 
relaxation in the energy transfer processes in photosyn- 
thetic antenna complexes has been studied theoretically 
in various variants [ 1 , H, B, H , HI, H E3, EH, El, EE Q 
ES ES El ES El M, MAM III ■ Essentially all works 
rely on the most commonly used approach which consid- 
ers the exciton-environment coupling as weak and treats 
it perturbatively [24| in order to derive some form of a 
time-local quantity for the time evolution of the system. 
This commonly results in or is equivalent to the famous 
Redfleld master equation El, EE El, El ■ Two central 
assumptions are made: (i) the system-environment cou- 
pling is sufficiently weak in order that a second-order 
perturbative description should be valid, and (ii) the 
bath-induced fluctuations of the excitonic energy lev- 
els are sufficiently fast so that they can be treated in 
a Markov sense. Only then, time-local evolution equa- 
tions result and the problem is sufficiently simple for fur- 
ther, often numerical, treatment. Put differently, it is 
assumed that the phonons which generate the fluctua- 
tions relax very quickly to their respective equilibrium 
state. The time scale thereby is given by the reorgani- 
zation energy Er which is assumed to be much larger 
than all other energy scales in the model. The reorga- 
nization energy is related to the cut-off frequency tu c of 
the bath spectral density J(w) and denotes the maximal 
frequency, up to which environmental modes exist. For 
the most important fluctuating electric dipole moments 
in the protein-solvent environment, the spectral distribu- 
tion of the bath modes has an Ohmic form @, m mi, i.e., 

J(u>) = 2awcxp(— uj/lu c ) with coupling constant a, which 
is related to the characteristic quantities of the solvent, 
such as the dipole moment, the frequency dependent di- 
electric constant and the Debye relaxation time [25[ . This 
implies Er = 2afku c . In order that the perturbative de- 
scription is valid, the intrinsic exciton dynamics (which is 
characterized by the energy gap A) has to be much slower 
than the bath reorganization processes, i.e, A <c u> c . 
Only then and with the additional weak coupling condi- 
tion Er <C u) c , all environmental influence can be treated 
to lowest order in the interaction and as effectively time 
local. 
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This, however, is typically not the case for the pho- 
tosynthetic excitation energy transfer. Typical energy 
scales in biomolecules are A ~ 0.2 — 2 meV for resonant 
energy transfer between red and green chromophores or 
A ~ 46 — 100 meV for the light-harvesting complex LH- 
II in the bacteria chlorophyll molecule in green sulphur 
bacteria [20l |. while the environmental energy scale are 
Huj c ~ 2 — 8 meV for water at room temperature. More- 
over, for typical free chromophores in water at room tem- 
perature, a ~ 0.1 — 1. This shows that one typically has 
to deal with cases where A < u> c , implying that any per- 
turbative approach to energy transfer in these structures 
becomes questionable. This has also been investigated in 
great detail by Ishizaki and Fleming [r| [HJ . They found 
that none of the variants of the perturbative approaches 
(Rcdficld equation in its full form, in the secular approx- 
imation and with the neglect of the imaginary part of the 
relaxation terms) can give a reliable picture [161 ] . In order 
to overcome the deficiencies, a reduced hierarchy equa- 
tion approach has been proposed [22j ■ At the same time 
numerical exact treatments were used [23] to describe the 
effect of slow environmental fluctuations. 

One of the experimentally best studied systems is the 
Fcnna-Matthews-Olson (FMO) protein (2fl [13] which 
serves in green sulfur bacteria such as Prosthecochloris 
aestuarii and Chlorobium tepidum as exciton conductor 
between the antenna complex and the reaction center. 
The FMO protein is a trimer where the three identical 
subunits consist of seven chromophoric sites which are 
fully characterized [28j. Typical transport times for ex- 
citons through FMO are of the order of picoseconds. Re- 
cent experiments [1^] on the FMO complex show that the 
exciton transport exhibits quantum coherent oscillations 
up to 660 femtoseconds. This raises the important ques- 
tion why the strong environmental fluctuations are in- 
sufficient to fully decohere the quantum transport. Sim- 
plifying the problem to a donor-acceptor model level it 
was shown [2^, [23[ that the slowness of the environment, 
namely lo c ~ A, sufficiently slows down decoherence for 
quantum oscillations to survive even for strong coupling 
to the environment. 

Having understood why quantum coherence survives 
rather strong couplings to environmental fluctuations 
raises the subsequent question whether photosynthesis 
benefits from quantum coherence either to improve ef- 
ficiency or to allow regulation. Detailed descriptions of 
transport through FMO, however, requires to go beyond 
a donor-acceptor model rendering the problem much 
more complicated. Therefore allmost all approaches to 
study this question rely on weak coupling approxima- 
tions JT3, M, E3, GJ, Q 01 GJ]- Ishizaki and Fleming 
[22I . l30| use a reduced hierarchy approach which they 
match for weak coupling to results from a Redfield ap- 
proach and for strong coupling to results from Forster 
transfer description. The actual parameter range of in- 
terest is in between. We show in the present paper that 
even at weak coupling perturbative weak coupling ap- 
proaches (like Redfield) break down and fail to describe 



the dynamics. 

We shortly review the Hamiltonian for a single chro- 
mophore in the next section |TT] before we once more in- 
troduce a donor-acceptor model in order to discuss the 
validity of perturbative approaches by comparing a low- 
est order perturbative approach with a numerical exact 
treatment. Using a resumed perturbative approach (RE- 
SPET) [H [S3, M, M, HI we analytically derive in sec- 
tion IIIII the time evolution of the population difference 
in our donor-acceptor model in lowest order in the sys- 
tem environment coupling. Thereby all non-Markovian 
corrections are determined in lowest order and they turn 
out to be negligible. In section [IV] the dynamics is solved 
using the numerical exa ct q uasi adiabatic propagator path 
integral (QUAPI) [H, H3 HI, HI E3- At weak cou- 
pling we find damped oscillatory behavior as expected 
from RESPET. However, for biomolecular environments 
at temperature T with A ~ lo c ~ T the numerically 
observed damping rate differs from the one-phonon rate 
of the analytical calculation by up to one order of mag- 
nitude. This difference increases with decreasing cut-off 
frequency ui c . This discrepancy cannot simply be ad- 
dressed to multi-phonon processes since the damping rate 
behaves strictly linear in the coupling and multi-phonon 
processes come only into play at stronger couplings. We 
conclude in section [V] accordingly that a weak coupling 
Markovian approximation qualitatively describes the dy- 
namics but fails quantitatively by up to an order of mag- 
nitude. This failure is nonperturbative in nature and 
can thus not be corrected by higher order perturbative 
approaches. Strictly spoken it invalidates all approaches 
to discuss the exciton dynamics in FMO which base on 
perturbative arguments. It even raises question marks on 
nonperturbative approaches which give at weak coupling 
comparable results to Redfield approaches. 



II. MODEL FOR BIOMOLECULAR ENERGY 
TRANSFER 

The simplest way to model a single chromophore (or 
pigment) is by describing it as a two-level system consist- 
ing of a ground and an excited state. When the electronic 
ground state is excited the electron is not free but local- 
ized by its attractive interaction with the hole it left. 
This dipolc electron-hole configuration forms an exciton. 
A formal description can be given as 

H = £y (1) 

where we introduced the Pauli matrix t z . Environmen- 
tal fluctuations will cause transitions between the ground 
and the excited state and will add a fluctuating energy. 
Experimentally it is known that the recombination time 
is of the order of nanoseconds whereas the complete en- 
ergy transfer through the complex is of the order of pi- 
coseconds. Thus the environmental fluctuations causing 
recombination are negligible. Describing the fluctuations 
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by harmonic oscillators, which couple linearly to the chro- 
mophore, results in the independent boson model for a 
single chromophore [1, H3, HH 



k 



(2) 



where we introduced the position and momentum oper- 
ators, qk and pk, of the mode with wave vector k and its 
coupling v k to the chromophore and fixed h = fee = 1 
which we keep below. 

We are not aiming at a description of the exciton trans- 
fer dynamics in the FMO complex but are interested in 
a thorough discussion of the validity of weak coupling 
perturbative approaches. For the sake of simplicity we 
thus restrict the model under consideration to just two 
chromophorcs. Since recombination is irrelevant for the 
transfer process of a single exciton we can restrict our 
investigation to two states \i) with the exciton at chro- 
mophore i. Furthermore neglecting the differing site en- 
ergies leads to a donor- acceptor Hamiltonian (5l[20l. l2lj 
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(3) 



with tunneling clement A between donor and acceptor 
and by introducing the Pauli matrices {I, a x , a y , a z } with 
a x = |1)(2| + |2)(1| and o z = |1)(1| - |2}(2| and 



A fc := fjfc(ri) - v k (r 2 ) . 



(4) 



Assuming a distance \r\ — r 2 \ > he/ A larger than the 
wave length of resonant modes with c being speed of 
sound the position dependence is negligible and the spec- 
trum 



k 



^ 2w fe 

k 



S(uk - to) ■ (5) 



The paramount approach to discuss the dynamics of a 
donor-acceptor model is the NIBA [4l|, [43]. It specifi- 
cally succeeds to describe the scaling limes, uj c — > 00, but 
is less trustable for bioenvironments with A ~ uj c . Be- 
sides numerical approaches [2^ the only other analytical 
rigorous approach is a weak coupling perturbative treat- 
ment which provides a treatment formally not limited to 
the scaling limes. 

We focus on the occupation difference between donor 
and acceptor P(t) = ((Jz)(t) with initially the exciton 
on the donor P(0) = 1 [4J|. Without coupling to the 
bath the occupation difference shows oscillations, P(t) = 
cos(At) with frequency A. The bath primarily causes 
damping of the oscillations. In order to determine the 
bath influence we calculate the time evolution operator 
U(t) with 



U(t) Pa = p{t) 
with the statistical operator p{t) obeying 
d tP {t) = -i[H,p(t)] =: Cp(t) 



(G) 



(7) 



with initial condition po and Liouvillc operator C in low- 
est relevant order in the coupling between spin and bath. 
We show in the next section that two additional weak 
contributions emerge whose residua are proportional to 
a and which are negligible under usual experimental con- 
ditions. 



III. ANALYTICAL APPROACH FOR WEAK 
SYSTEM-BATH COUPLING - RESPET 

For a perturbative approach in the coupling between 
spin and bath we split the Hamiltonian 

H = H -H SB with H SB = -^^A fc g fc (8) 



and likewise the Liouvillian C = Co + Csb- The time 
evolution operator obeys a Dyson equation 



U(t,t ) = tto(Mo) + / dsUo(t,t )£ SB U (s,t ) 

Jto 

ft ps 

ds / ds'Uo(t,t )£sBUo(s,s')£sBU(s / \t ) (9) 



with Uo(tito) the time evolution of the uncoupled spin 
and bath. We are interested in the dynamics of the spin 
alone allowing us to integrate out the bath degrees of 
freedom in the Dyson equation resulting in 

Ues(t,t ) = U S {t,t ) 



ds / ds'U s (t,t )M(s,s')U eS {s',t )(10) 

to Jto 



with the memory kernel M(s,s') = (LsbUq{s,s')Csb)b 
where (-)b = Tr{-ps} with the initial bath statisti- 
cal operator ps assuming a factorizablc initial condition 
Po = Ps ® Pb- Herein it is essential that the linear order 
term (Csb)b = vanishes. The memory kernel is given 
in lowest order (1-loop approximation) in the coupling 
between system and bath only 44 1 . 

With the Laplace transformation defined as 



/(*) 



i I dte lzt f(t) 
'0 



and 



1 r 

m " Si/., 



7M 



the Dyson equation translates into 

U cS {z) = (Uj\z) +M(z))~ 1 



(11) 



(12) 



(13) 



For a single two level system or spin the space for 
the time evolution operator U e g is spanned by basis 
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{11, a x , <7 y , a z } and the operators U c g as well as Us and 
M. can be expressed as 4®4 matrices 



U s {z) 



-z 






-z~ 







z 2 -A 2 

-iA 
z 2 -A 2 



iA 



(14) 



-A? 









\ 


-(j){z) -ip(z) 







-2T(z) 










o / 



M(z) 



with 

<j>{z) = B R {z - A) - B R (z + A) 

V>(z) = 2{B c (z- A)+B c (z + A)} 

T(z) = 2B c (z) 

and the bath correlation function 



(15) 



Bc(z) 



dujJ(cu)coth(/3u>/2)- 



and the bath response function 



Br(z) 



dhj2J(uj)- 



(16) 



(17) 



with J(uj) = J2k(^l/ 2rnuj k)S(uJ - LOk) = 
2awexp(— cj/w c ). 

Finally we obtain the effective time evolution operator 



U c fi{z) 



( zl ° \ 

a(z) b(z) 

z -iA 

n p( z ) p( z ) 

U _iA_ z+2T(z) 

viz) v(z) I 



(18) 



with 



a(z) 
b(z) 



-{<P{z)M>{z)){{z + iP{z))- 1 - z~ 1 } 

-(z+nzjr 1 



and 

p(z) = T 2 (z) - {z + A + T(z)}{z - A + T(z)} (19) 

We focus on the occupation difference of the two level 
system 



P(t) := (a z (t)) s 



(20) 



where (-)s = Trj-pg} with the systems initial statistical 
operator ps = \ (11 + cr 2 ) defining P(0) = 1 as initial 
condition. With the definition of correlation operators <7j 
by (fjO-j = 7}{(Ti(jj +aj<Ji} one can express the occupation 
difference by 



P(z) = (a z U oS (z))< 



z + 2T{z) 
p(z) 



(21) 



Thus we have formally solved the dynamics of the two 
level system. What is left is the Laplace back transfor- 
mation: 



P(t) = ^~ 



dze 



: z + 2T(z) 
p(z) 



(22) 



The integration is closed in the lower half plane and thus 
we need to know the branch cuts and poles of P(z) in the 
lower half plane. Keep in mind that the Laplace trans- 
forms are analytical in the upper half plane by definition. 

The integral representation of the bath correlation (fTB|) 
is not valid in the lower half plane since it has a branch 
cut along the real axis. In order to continue the function 
into the lower plane we split it according to Be = Bq + 
B^, with 



du>J(uj)coth{/3uj/2) (23) 

Z^fUJ 



where B^ has its branch cut only from to ±oo and thus 
can trivially be continued for Re(z) ^ 0. For Re(z) > 
we need to shift the path of integration in the integral 
representation of Bq in order to continue the function 
into the lower half plane leading to 



OO ^ 

dtoJ(uj) coth(/3w/2) 

o z T w 



+i\nJ{±z) coth(/3(±z)/2)9 (Im(z) ^ 0) (24) 

with the Heavyside function 6(-). Since the integrand, 
i.e. the hyperbolic cotangent, has poles on the negative 
imaginary axis the path of the integrand can maximally 
be shifted to run along the negative imaginary axis. Thus 
the bath correlation function can be presented having a 
branch cut along the negative imaginary axis and as well 
poles on it due to the term coth(/3(±z)/2). 



A. System poles 

With this prerequisite we can tackle the Laplace back 
transformation. We expect generally three contributions. 
First without coupling to the bath the dynamics is gov- 
erned by the two system poles z± = ±A leading to 
P(t) = cos(At). At weak coupling we expect the main 
behavior still governed by these two poles but we further 
expect corrections to the frequency, a decay rate and the 
prefactor diminishing this contribution to the total dy- 
namics since the branch cut as well as the poles of the 
bath correlation function will result in additional contri- 
bution. 

Without going into details the system pole contribu- 
tion becomes 

P s (t) = + 3Jm(r(A)) + Rc ^ A )^ cos(At)e- r " t 
d z Re(r(A)) + Im( ^ (A)) ) sm(A;)e- r e25) 



with A = A + Re(r(A)) and T = Im(T(A)) = 
jirJ(A) coth(/3A/2). A standard Markov approximation 
stops at this point since it focuses on the dynamics result- 
ing from the system poles. Within the presented weak 
coupling approach the other contributions can also be 
determined. 



B. Bath poles 

The second contribution comes from the poles of the 
bath correlation function since around them T(z) will 
become very large and accordingly 



p(z) 



z 



2T(z)z 



(26) 



will have zeros. The poles are at — it n = —i2imT with 
?i > an integer. Focusing first on Re(z) > we can 
safely assume that the integral part in the continuated 
bath correlation function will be small at the pole of the 
second part and thus we obtain 



-it r . 



|7riJ(Im(r(A))) 



2T 



(27) 



which leads us to zeros of p(z n ) at z n = —it n + e+ with 
+ _ irt n J(-it n )2T 



tl + A^ 



(28) 



Depending on the ratio t n /cu c the real part of e+ turns 
out to be positive or negative. Whenever it is negative 
for a certain n this pole does not contribute since we 
restricted the investigation to Re(z) > 0. Similarly we 
get 



nt n J(it n )2T 
€ + A2 



(29) 



which only results in a pole when Rc(e~) < 0. These con- 
ditions arc only fulfilled for both identically ensuring that 
all poles come in pairs resulting in a real time-dependent 
function. The residua of these contributions can be esti- 
mated by 



Res 



fz + 2T(z)} 




= Re ^ 


I P(z) J 





A 2 T 



i\j{-it n ) I = 



finally leading us to 



P B {t) = ^%„cos(Rc(e+)i)e~ 



(t„+Im(e+)t 



(30) 
(31) 



With increasing n the weights v n decrease and the re- 
laxation times increase for these terms. In the case 
T ^ A ~ oj Ci we are interested in, only the very first 
few n's are of any relevance. The relaxation time is of 
the same order as the tunneling element A and thus this 
term is fast decaying and weak due to v n to begin with. 
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FIG. 1: Fhe occupation difference P(t) versus time for T = 
uj c — 1.5A. Symbols represent data (most points not printed) 
obtained by QUAPI and the full lines are fits to the data. 



Its frequency, however, is small too, thus when T ~ A 
on a time scale of A -1 this contribution will be present 
as an ovcrdamped behavior for not too weak couplings. 
This contribution is clearly a bath generated dynamics 
which is non-Markovian in nature. 



C. Branch cut 

The third contribution to the occupation difference 
P(t) stems from the branch cut of the bath correlation 
function which generates a similar branch cut along the 
negative imaginary axis for P(z). Firstly we have to cal- 
culate the jump at the branch cut 



S(z 



—is) = lim{P(?7 — is) — P(— T] 
ri— >o 



iS )} (32) 



Close to the poles T(z) dominates and the jump van- 
ishes. As long as we are far from the poles of T(z) we 
can approximate 



S(z = —is) ~ — 2-7T 



Re{J(-is)}cot(/3s/2) 



and then obtain 



P G (t) 



-1 
2^ 



dser st S(-is) 



(33) 



(34) 



We are interested in times t > A -1 and T ~ A ~ uj c 
and thus the exponential in the integral allows us to re- 
strict ourselves to s < A and to approximate S(—is) ~ 
-47ras(2T/A 2 /w c ) and thus 



a T 1 
4a 

to c (At) 2 



(35) 



This contribution could actually be rather large when 
(T/co c ) > 1 and its decay is algebraic and thus slow com- 
pared to the exponential decay of the other terms. How- 
ever, we are restricted to times t ^> A^ 1 rendering its 
contribution small. 
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FIG. 2: The decoherence rate r num versus coupling constant 
a at T = A = u) c including error bars for the numerical 
procedure. The line is a linear fit to the data. 



All in all, we have three contributions to the dynamics 
of the two level system. We denote them in the following 
as the system pole, the bath pole and the bath branch 
cut dynamics. The contribution of the later two is pro- 
portional to the coupling constant a and thus small in 
weak coupling. These are the non-Markovian parts of the 
dynamic. An estimate of the three parts for A = T = uj c 
at a time t\ = 2 A -1 is obtained as 



a 



E 



Trnae 



-4nn 



nae 



-4tt 



< 2(47r 2 n 2 + 1) 2(4tt 2 + 1) 



P s (h) < exp(— roe -1 coth(i)) ~ e 



-2.5a 



and shows that the non Markovian parts arc only small 
corrections at weak coupling. Thus, one would conclude 
that at weak coupling Markovian dynamics and a one- 
phonon process dominate. In the next section we deter- 
mine the occupation difference using the numerical exact 
QUAPI [H, [33, [H, HI, 53] approach in order to test the 
validity range of the weak coupling approach. 



IV. NUMERICAL APPROACH: 
RESULTS 



EXACT 



Usin g th e q uasi adiabatic path integral (QUAPI) [36l . 
l37l. I3RI39I l4fjj|. a numerically exact approach, we deter- 
mined the time dependent occupation difference P{t) = 
(cr z )(i) of the donor- acceptor system coupled to envi- 
ronmental fluctuations with Ohmic spectrum and A ~ 
fceT ~ w c - The occupation difference P(t) shows weakly 
damped oscillations which we fitted with 



P(t) = {cos(A nmn £)+Msin(A 

num 

t)}e 



(36) 



using the tunneling frequency A num , the damping or de- 
coherence rate r num and a prefactor \i as fitting param- 
eters. 
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FIG. 3: The ratio r n um/ri_ p hon versus coupling constant a 



According to our weak coupling analytical result the 
prefactor [i is proportional to the coupling a. The fitting 
accuracy, however, is rather insensitive to [i. Besides it 
is small for small a and for large a the occupation differ- 
ence is dominated by the exponential decay. Therefore we 
fixed n = r num /A num as known for large iv c [HI, [HI. The 
modification of the tunneling frequency, |A num — A|/A, 
due to the coupling to the environmental fluctuations is 
less than 5% in the whole studied parameter range. This 
illustrates that for oj c < A the effect of the fluctuations 
on A cannot be included by a mere parameter renormal- 
ization A — » A e ff [4l|, [4^] , as it is the case in the regime 
lu c 3> A. The main relevant effect of the environment is 
the emergence of a decoherence rate which we discuss in 
the following. 

The chosen fit function was sufficient to describe all 
data discussed. Fig. [T] shows exemplary two data sets 
with T = uj c = 1.5A for a weak coupling a = 0.1/(27r) 
and a medium coupling a = 1.5 /(2n) with the respective 
fits. 

Now we focus solely on a discussion of the decoher- 
ence rate. Fig. [5] plots the decoherence rate at tem- 
perature T = A and for a cut-off frequency u c = A 
versus coupling a. The line in fig. [5] is a fit to the 
numerical data giving clear evidence for the linear de- 
pendence of the decoherence rate on coupling strength: 
r nnm = 0.246(27r)aA. Our weak coupling result 
predicted T = \ cxp(-A/cj c ) coth(A/2fc B T)27raA ~ 
0.199(27r)aA. Thus the numerically extracted decoher- 
ence rate is about 25% larger than the weak coupling ex- 
pression predicts for all investigated coupling strengths. 

At first glance one might attribute these differences to 
multi-phonon processes which are beyond the analytical 
onc-phonon rate. However, multi-phonon processes gen- 
erally are not linear in the coupling a and accordingly 
the decoherence rate should not be linear in a. In or- 
der to see the onset of multi-phonon processes at higher 
couplings Fig. [3] shows the ratio of the numerically ob- 
served decoherence rate to the 1-phonon rate versus the 
coupling strength a for three temperatures and cut-off 



7 



= 10 



a 


1 1 

□ 

: O 


i i i i ii | 


i i i i i i ii | 

O T = 0.2 A " 
□ T = 0.5 A 
O T = 5 A _ 








- 


s 

3 








c 

u 








1 


i i 


o g 

1 


i 



0.1 1 10 



cut-off frequency © / A 

FIG. 4: The ratio r n um/ri_ p hon versus cut-off frequency uj c 
for three different temperatures for a fixed weak coupling a = 
0.1/ (27r). 



frequencies. At weak coupling, a < 0.2/(2%) the ratio 
is constant. At medium couplings 0.2/(271") < a < 0.5 
the ratio slightly increases with coupling showing the on- 
set of multi phonon processes (higher order processes in 
a) [45[. Most importantly, the deviations of the ratio 
from 1 arc independent of the coupling strength for weak 
couplings. The deviations increase slightly with increas- 
ing temperature and rather strongly with decreasing cut- 
off frequency u> c . 

Fig. [2] and Fig. [3] together show that in the parame- 
ter regime A ~ k^T ~ uj c relevant for bioenvironments 
the weak coupling approach fails even at lowest coupling. 
Since the deviations between the 1-phonon rate and the 
numerical exact one is independent of the coupling a (for 
weak coupling) it is not a result of higher order processes 
and thus perturbatively not tractable. Accordingly all 
perturbative approaches (even in higher order) are bound 
to fail to describe the dynamics correctly beyond a mere 
qualitative description. 

In order to show the full extent of this failure Fig. [4] 
plots the ratio r num /ri_ p hon versus cut-off frequency lo c 
for three different temperatures for a fixed weak coupling 
a = 0.1/(2%). For large cut-off frequencies oj c 3> A the 
ratio tends to 1 and the onc-phonon process accurately 
describes decohering dynamics. With decreasing cut-off 
frequency the ratio increases strongly and at lo c = A/5 
the numerical exact rate is an order of magnitude larger 
than the one-phonon process predicts. 

Beyond our RESPET result we compare the numeri- 
cal results with, one might be tempted to use available 
more sophisticated analytical treatments [4l|, [42| . These 
give in the limit of weak coupling expressions for the de- 
coherence rate Tq formally identical to our result (and 
thus Fermis golden rule) but the tunneling element A is 
exchanged by a renormalized one, A e fi oc (A/uj c ) a ^ 1 ~ a \ 
depending on the cut-off frequency. At first we should 
remark that this renormalization decreases the decohcr- 
ence rate for decreasing lo c (in comparisson to ri_ p i lon ) 



which is in clear contrast to our finding and secondly 
the adiabatic renormalization leading to the expression 
of A e g is justified only in the limit Acu c whereas our 
observed deviations occur in the opposite limit. 



V. CONCLUSION 

Exciton transport in photosynthesis is strongly dis- 
turbed by the fluctuations of the dipolar solvent for 
the biomolecules building the light-harvesting complexes. 
However, the fact that these fluctuations are slow, uj c ~ 
A, results in far longer decoherence times than known 
for solid state excitons where lo c 3> A and finally allows 
for enhanced quantum coherence in such systems [23| . 
To study the extent and relevance of quantum coher- 
ence on exciton transport, however, systems like the 
FMO complex must be studied beyond a donor-acceptor 
model. Numerically exact treatments, like QUAPI or 
QMC, become rather time consuming when including all 
7 chromophores and one is restricted to analytical ap- 
proximations. Slow but strongly coupled environmental 
fluctuations cause strongly non-Markovian dynamics and 
thus approximative analytical approaches must be able 
to include such effects. To verify such methods they are 
typically matched to weak coupling perturbative results 
and to Forster transfer dynamics in the weak tunneling 
limit EH]. 

We have determined the population difference in a 
donor-acceptor model in lowest order in the coupling 
between donor-acceptor and environment including all 
non-Markovian contributions and compare it to results 
from a numerical exact treatment using QUAPI. We show 
that in the experimentally relevant parameter regime, 
A ~ k^T ~ uj c , where there is no small parameter, the 
decoherence rate strongly differs from the one-phonon 
rate and the deviation increases with decreasing cut-off 
frequency u c . This was shown for a weak coupling where 
the decoherence rate varies strictly linear with the cou- 
pling. Multi-phonon contributions showing non-linear 
dependence on the coupling arc found for larger couplings 
as well. The deviations also increase with increasing tem- 
perature. 

Thus we conclude that analytical approximations to 
discuss the influence of bioenvironmental slow fluctua- 
tions on quantum systems cannot be verified in the weak 
coupling limit by comparing to weak coupling perturba- 
tive results. 

The increasing deviations with decreasing cut-off fre- 
quency lj c points to strong nonperturbative dynamics for 
small cut-off frequencies. We interpret this behavior as 
the breakdown of a picture where the slow environment 
acts as a bath which only exchanges energy with the 
quantum system. The fact, however, that the dynam- 
ics can still be described by damped coherent oscillations 
with a decoherence rate r num varying linear with cou- 
pling constant a shows that at weak coupling a slow bath 
renormalizes itself to an effective bath whose influence 
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can be described by an effective one-phonon process with 
renormalized bath parameters. An analytical description 
of this renormalized bath would be needed to base ap- 
proximations on which allow to study the exciton trans- 



port dynamics in photosynthctic light-harvesting com- 
plexes. 

We gratefully acknowledge support by the Excellence 
Initiative of the German Federal and State Governments. 
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